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Precision measurements in charm pliysics offer a window into a unique sector of potential New Physics interac- 
tions. LHCb is well equipped to take advantage of the enormous production cross-section of charm mesons in pp 
collisions at y/s = 7 TeV. The measurement of the D'^ — mixing parameters and the search for CP-violation 
in the charm sector are key physics goals of the LHCb programme. The first CP violation measurements in the 
charm sector, with 37 pb~^ of data collected in 2010, are discussed. The study of — > K~ decays 
shows no indication of CP violation. The measurement of the proper time asymmetry in the time dependent 
analysis of D° K' K+ and D° -5- K' K+ is evaluated to be Ar = (-5.9 ± 5.9stat ± 2.1sy3t)- The difference 
of CP asymmetry in the time integrated rates of D*^ — > K~ and D^^ — > 7r~7r+ decays is measured to be 
(-0.28 ± O.lOstat ± 0.25sast)%- 



1. Introduction 

LHCb [T], an experiment at the Large Hadron Colhder (LHC), is dedicated to the study of 6 and c flavour 
physics. The abundance of charm particles produced in LHC offers an unprecedented opportunity for high 
precision measurements in the charm sector, including measurements of CP violation and — mixing. The 
high performance of LHCb detectors allows this potential to be fully exploited. 

The detector is a single-arm forward spectrometer covering the geometrical region where heavy flavour par- 
ticles, at LHC energy, are mostly produced. A silicon micro-strip vertex detector (VELO) provides, with high 
precision, the position of the primary vertex and those of the decay of long-lived particles. Other elements of 
the LHCb tracking system include a silicon strip detector (TT) located in front of a dipole magnet and three 
station detector downstream of the magnet, composed of a silicon micro-strip detector (IT) in the inner part 
and by straw drift chambers (OT) in the outer region. Charged hadron identification is made through two 
ring-imaging Cherenkov detectors (RICH). The calorimeter system identifies high transverse energy hadron, 
electron and photon candidates and provides information for the trigger. The particle identification system is 
completed by five muon stations that provide fast information for the trigger and muon tagging. 

CP violation in D decay processes has not yet been observed. In the SM, indirect CP violation in the charm 
sector is expected to be highly suppressed, less than 0(10"^), and universal between CP eigenstates. While, 
direct CP violation can be larger in SM dependent on the final state: CKM dynamics can produce direct CPV 
asymmetries in single Cabibbo suppressed decays of the order of 10~'^ or less [2]. Both asymmetries can be 
enhanced by New Physics in principle up to 0(1%) [3]. 

In 2010 LHCb recorded a total integrated luminosity of 37 pb^^. This provides a charm sample large enough 
to be able already to make several competitive measurements. The expected integrated luminosity of more 
than 1 fb~^ foreseen in 2011 will offer the opportunity to improve the world knowledge of D mixing and CP 
violation. The results of search for direct and indirect CP violation in the charm sector on data taken in 2010 
by LHCb are presented here. In particular, the search for direct CP violation in singly Cabibbo suppressed 
(SCS) decay D'^ — ^ K~ K^t:^ , the measurement of indirect CP violation in mixing in two body hadronic 
charm decays, and the search for CP asymmetry in the time integrated rates of D mesons into 2 body SCS 
decays are illustrated. 



2. Search for CP violation \x\ D"^ ^ K K^-n^ decays 

An independent analysis is performed to search for direct CP violation in the singly Cabibbo suppressed 
decay D'^ — K~ K^-k'^ . The search consists of a direct comparison between the and the D~ Dalitz plots 
on a bin-by-bin basis. The Dalitz plot is divided into bins and for each bin a local CP asymmetry variable is 
defined: 

Srp ^ , with a = ) -, (1) 

^ ^m{D+)+am'{D-) Ntot{D-y ^' 

where N^{D^) and N'^{D~) are the numbers of decays in the ith bin and a is the ratio between the total 
and yields. The parameter a is a correction to account for global asymmetries that are constant across 
the Dalitz plot. 
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Figure 1: Mass spectra of K K^tv'^ . The signal mass windows and sidebands (lower, middle and upper) are labelled. 



In the absence of local asymmetries, the S^p values are distributed according to a Gaussian distribution 
with zero mean and unit width. CPV signals are, therefore, deviations from this behaviour. The comparison 
between the Z?"*" and the D~ Dalitz plots is made by a test [4]. The is defined as = J2{^cp)^ ^^'^ the 
number of degrees of freedom (ndof) is the number of bins minus one. Hence the probability value (p-value) 
measures the confidence level that the difference between Z)"*" and D~ Dalitz plots is driven only by statistical 
fluctuations. 

Different binning schemes are considered to obtain the highest sensitivity to various types of CPV. The bin 
scheme was optimized taking into account that we have no sensitivity if CP asymmetries change sign within a 
bin and we have a reduced sensitivity if only a small part of a large bin has any CP violation in it. 

The technique relies on careful accounting for local asymmetries that could be induced by sources such as 
the different production mechanisms for D+ and D~ , the difference in the K-nucleon inelastic cross-section, 
differences in the reconstruction or trigger efficiencies, left-right detector asymmetries, etc. The existence of 
these local asymmetries are investigated using the Cabibbo favoured control channels, D"*" K~'K'^n~^ and 
_D+ — > K^TT^ . No CP violation is expected in these channels. The first control mode, _D+ K^tt^tt^, 
has an order of magnitude greater branching ratio than the Cabibbo suppressed signal mode and is more 
sensitive to detector effects since there is no cancellation between and K~ . The second control mode, 
Df — > K^TT'^ , is similar to the signal mode in terms of resonant structure, statistics, kinematics, detector 
effects and backgrounds. Similarly, the method is also applied in the sidebands (shown in Fig. [I]) of the second 
control channel to investigate possible asymmetries due to the contamination of the background. Another source 
of asymmetries could come from a charge asymmetry from the parent B in the B — >■ D{K~ K~^Tr^)X decays. 
The effect of secondary charm is investigated by dividing the data set by the impact parameter-'^ (IP) significance 
ixjp) to have samples with different contamination from secondary charm. All these tests are fully consistent 
with no asymmetry, thus the method is determined to be very robust against systematic effects. 

The data sample used in this analysis corresponds to approximately 35 pb"""^ collected in 20f0. The signal 



Table I: The p-values for consistency with no CPV for the D° — >■ K K'^n'^ decay mode for data with different magnet 
polarities 



Magnet Polarity 


p-value 


Up 


6.0% 


Down 


28.5% 


Combined 


12.7% 



^ The IP is the minimum distance of approach with respect to the primary vertex. The X^p is formed by using the hypothesis 
that the IP is equal to zero. 
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Figure 2: Distribution of S^p in the Dalitz plot for the Uniform scheme. 



sample consist of about 370k candidates. The global asymmetries parameter (a) is measured to be 0.984±0.003. 
Fig. [l] shows the invariant mass of K~ K^tt^ for the analysed data sample. The following two bin schemes are 
used on the signal data: the first scheme (Uniform) uses an uniform grid of equal size bins; the second type 
(Adaptive) takes into account the non uniform event distribution due to the 7r+ and i?*(892)°i<'~ modes. 
This second scheme has bins of variable size, aiming for a uniform population in all the bins. For each bin the 
significance S}jp of the difference in 2?+ and population is computed as defined in Eq. [ijand the x^/ndo/ 
is calculated to obtain the p-value. The data with opposite magnet polarities are combined to cancel left-right 
asymmetries. 

The obtained p- values, summarized in Table|lj indicate no evidence for CPV [5 . This result is also supported 
by the result of the Gaussian fit of S^^p that have mean and width consistent with and 1, respectively. The 
results for the Uniform bin scheme are shown in Fig. [2] 



3. Measurement of indirect CP violation in mixing 

A measurement of the indirect CP violation in mixing can be performed in the study of two-body hadronic 
charm decays. It can be evaluated by the asymmetry of the proper-time (r) of flavour-tagged decays: 

T (Z)o ^ K-K+) -t(D" ^ K-K+) 1 ^ 

Ar = — \ ) — r = - Am V cosd) — X sin0, (2) 

T {D^ ^ K- K+) + T {D^ ^ K- K+) 2 ^ ^ ^' 

where x = and y — ^ are the mixing parameters and cj) is the CP violating weak phase. Am is defined by 
the parameterization —\ q/p |^^= 1 ± Am with the assumption that R^ is close to unity and where q and 
p are parameters that define the mass eigenstates in terms of the flavour eigenstates. 

A measurement of Ar differing significantly from zero would be a measurement of indirect CP violation as it 
requires a non-zero value for Am or (j). 

The signal yield and the background contribution are extracted from fits to the reconstructed invariant mass 
alone. Due to the abundance of charm decays, the selection has been designed to achieve maximal purity, with 
a background rate of the order of a few percent. The main component of the background is due to the secondary 
charm, i.e. D mesons produced from b hadron decays. This kind of background is not distinguishable by the 
invariant mass distribution. The secondaries have larger impact parameter with respect to the primary vertex 
than the prompts as a secondary D no longer has to point back to the primary vertex. Thus this background 
can be reduced by a selection based on the topology but it can not be completely suppressed. Hence a statistical 
separation is required. We use the variable In {xjp) i because it is an easier quantity to parameterise than IP 
directly. For the secondary charm, the IP depends on the B flight distance. The form of the consequent 
In {xjp) dependence on proper-time is extracted from the simulation, and the parameters of this dependence 
are evaluated in the fit procedure. 

Since this analysis is sensitive to the proper time dependence of the acceptance, particular attention is paid to 
requirements that could bias this distribution. A correction of these lifetime biasing effects is needed to properly 
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Figure 3: In (x/p) fit projection of D" — >■ K~ K'^ candidates. The data are shown as points, the total fit (blue), the 
prompt signal (red), and the secondary signal (pink). 

extract via absolute lifetime measurements. The heavy flavour selection implies some criteria which bias the 
measured proper time distribution. These biasing selections are unavoidable and have to be applied already at 
trigger level to suppress background from the large number of particles produced promptly in the proton-proton 
collisions. One example in the LHCb selection of heavy flavour decays is the request of a minimum impact 
parameter of the daughters. A correction of these lifetime biases has to be computed to properly evaluate the 
lifetime to determine Ay- 

This analysis uses a data driven approach to evaluate the proper time acceptance, that describes the selection 
efficiency as a function of the proper time. The method evaluates the proper-time acceptance on a per 
event basis by a so-called 'swimming' algorithm, which was originally developed at CDF [S] [7| and is now being 
applied at LHCb [8 lOj. fn this method the acceptance function is evaluated by moving the primary vertices 
(and thus varying the proper time for the candidate). For each primary vertex position the software trigger 
decision and the offline selection are re-evaluated. Consequently the proper-time acceptance function for each 
event is determined as a sum of step functions, indicating when the event would be selected or not selected. 

The measurement of is performed via absolute lifetime measurements obtained by a simultaneous fit of 
proper time and In (x/p) including the acceptance function evaluated by the swimming method. Fig. |3] shows 
an example of the projection of the In(xjp) for _D" K^K^ decays. 

The measurement is based on a data sample equivalent to 28 ± 3 pb~^ of data taken in 2010. The number 
of candidates selected is about 15k for each flavour tag, and D'^. The flavour tagging of D'^ decays is done 
by reconstructing the decay D*^ — > D'^tt^, where the charge of the slow pion (tTs) determines the flavour of the 
-D" at production. 

The method was validated on a control measurement using decays to the Cabibbo favoured decay D — Kn. 
The combinatorial background contribution is ^ 1 % for the control channel and ^ 3 % for — > KK decays. 
In the 2010 data sample only low statistics were available in the mass sidebands. Hence it was not possible to 
model the background shape, and the background contribution is neglected in the time dependent fit and it is 
taken into account in the systematic uncertainties. 

In the control channel the result for the lifetimes, averaged between 13° — )■ K~7t~^ and _D'^ — K~^Tr~, is 
t(D°) — 410.3 ± 0.9 fs, where the uncertainty is statistical only. This is in agreement with the current world 
average [T2]. The lifetime asymmetry has been determined as Ap'^ = (—0.9 ± 2.2stat ± ^-Qsyst) ■ 10~^, which is 
consistent with zero in accordance with the expectation. 

The measured lifetime is an effective lifetime since the fitted distribution includes also mistagged events, in 
which the is associated with a random slow pion. The mistag rates are assumed to be independent of the 
final state and are extracted from the favoured D Ktt decays which offer higher statistical precision. This 
rate is evaluated by the fit of the difference between the mass of D* and (Am) to be 1.8%. This has been 
neglected in the control channel as it is very small, but it is applied in the evaluation of Ar- The distribution 
of the Am between the D* and D° is shown in Fig. |4j 

The results of the lifetime flt of — ?► K^K^ and K^K^ are shown in Fig. [s] The asymmetry is 

evaluated from these lifetimes to be |TT]: 
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Ar = (-5.9 ± 5.9stat ± ^.Isyst) ■ 10"^ 



(3) 
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Figure 4: Am mass difference between the reconstructed D* and D candidates. The data are shown as points, the total 
fit as a solid line and the random background as dashed line. 




Figure 5: Lifetime fit projection of D^^ — >■ K candidates on the left and D° — > K candidates on the right on a 
logarithmic scale. The data are shown as points, the total fit (blue), the prompt signal (red), and the secondary signal 
(pink). 



This result is consistent with zero and hence shows no evidence of CP violation and is in agreement with 
the current world average |12) . The main contributions to the systematic error are due to neglecting the 
combinatorial background and to the separation of prompt and secondary charm decays. The systematic 
uncertainty is expected to be significantly reduced by an improved treatment of the background events, which 
will be possible for the data taken in 2011. 



4. Search for CP asymmetry in the time integrated decay rates of D mesons 

LHCb is searching for evidence of new sources of CP asymmetry in the time-integrated decay rates of D 
mesons. The asymmetry is defined as 



^RAw,fy ^ N{D^^f)-N (D° ^ /(/) ^ N {D*+ D^{f)n+) - N {P*- ^ D^{f)n-) 
N{DO ^ f) + N{D° ^ fif) N{D*+ ^ D°{f)n+)+N{D*- ^ DO{f)Tr~) 



where N{X) refers to the number of reconstructed events of decay X after background subtraction. The raw 
time integrated asymmetries of and Z?" decays are considered separately using the D* decay slow pion 
tagging method explained above. 

The raw asymmetries may be written as a sum of various components, coming from both physics and detector 



6 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



2200r 

\ 2000 

« i8oo=-^"'ri^i„,,„ 

lo ^gQQ ^ Preliminary 

iS! 14001-^='"^"°^"' 

~ 12001- 
B 1000 1- 

400 r 
200 H 




2.006 2.008 2.01 2.012 2.014 2.016 2.018 
m(KK™)-m(KK)+m(D\ „^ (GeV) 




2.006 2.008 2.01 2.012 2.014 2.016 2.018 



m(n;iui)-m(mi)+m(D (GeV) 



2.006 2.008 2.01 2.012 2.014 2.016 2.018 



2.006 2.008 2.01 2.012 2.014 2.016 2.018 



Figure 6: Example fits to the mass difFerence spectra for tagged candidates for a subset of the data. The decay is 
reconstructed in the final states K~ K'^ (left) and ■K~n'^ (right). The normalized residuals for each bin are shown below. 



effects: 

A^f^ifY = Acp{f) + Anif) + Ap{D*), (5) 

where Acp{f) is the physics CP asymmetry, Ani^f) the detection asymmetry of the Ao{t:s) the detection 
asymmetry of the soft pion and Ap{D*) the production asymmetry. 

Taking the asymmetry difference of the two final state {ApAwifY ~ -^RAwU'Y) the production and soft 
pion detection asymmetries will cancel. Moreover, for a two body decay of a spin-0 particle to a self-conjugate 
final state, there is no D° detector efficiency asymmetry contribution, i.e. A]j(K^ K^) = A]j{tt^tt^) = 0. 
Due to possible production and detection asymmetries, the measurement of time-integrated CP asymmetry 
independently in D'^ -> KK and D'^ — ^ tttt is challenging. 

We can however measure the difference in time-integrated CP asymmetry AAcp between — > KK and 

D" TTTT. 

A^f^{KKY = Acp{KK) + Ad{^,) + Ap{D*+), (6) 

Ag^^(7r7r)* = Acpi^ii) + AoiiTs) + Ap{D*+), (7) 

Araw{KKY - A^f^i^^Y = Acp{KK) + Acp(w). (8) 

No dependence remains on production or detection efficiencies, so this observable is extremely robust against 
systematic biases. 

In a proton-proton collider machine the production of heavy-flavour hadrons need not be CP symmetric in 
a given region of phase space. Possible variations of both selection efficiency and production and detection 
asymmetry as a function of pr and rj could generate second-order yield asymmetries that do not cancel out 
m our formahsmjs] A^f^ extraction is performed in bins of rj and pt chosen such that the statistics are 
approximately constant within each bin. The binning is chosen to take into account the potential variation of 
production or detection asymmetries in these variables that differ for the two final states, such as those that 
may be induced in the selection by e.g. particle identification requirements. 

A binned maximum likelihood fit to the spectrum of the mass difference between D* and D'^ is used to 
evaluate the yields. Examples of the fit are shown in FiglH] The data sample has an integrate luminosity of 
37 pb~^. The total signal yield is 116k tagged Z?" — K~K~ and 36k tagged — )■ 7r~7r+ . The background of 
mis-reconstructed D* decays that peaks in the mass difference is estimated from the mass sideband to be at the 
sub-percent level. The effect enters the asymmetry calculation at second order, 0(10^"'), and can be neglected. 

Systematic uncertainties are assigned by repeating the analysis with an alternative description of the mass 
spectra lineshapes; with different fit windows for the mass; with all candidates, choosing one candidate 
randomly in events containing multiple candidates; and comparing with the result obtained with no (ptiV) 
binning. The full change in result is taken as a systematic uncertainty and all uncertainties are added in 
quadrature. No source of limiting systematic bias has been identified. These uncertainties are expected to be 
reduced by exploiting the much larger statistics that will be available. 

A value of AAqp is determined in each measurement bin using the result from Ap^wi^^^ ^'^^Y ^^'^ 
ApAwiT^^T^^Y ■ These values are found to be consistent throughout the {pt,ii) space, as well as for the 
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two trigger periods and for both settings of the magnet polarity. A weighted average is therefore performed 
to yield the result A^cp = (-0.28 ± 0.70 ± 0.25)% [13 . This result is approaching the sensitivity of CPV 
measurements performed by the B-factories in these decay modes[Tll [TS], but not yet at the level of CDF's 
recent measurement |I6j. 

The time- integrated CP asymmetry A^cp between the final states K^K^ and D'^ — 7r~7r+ has two 

contributions: a direct and an indirect component, a^p and a^p, respectively. The indirect component may 
be assumed to be the same for both final states as it originates in the common box diagram. However, its time 
dependence has to be taken into account, leading to a non-cancellation if the two final states are reconstructed 
with a different mean proper time. While the direct component is different, in general, for different final state. 
Thus, the physics asymmetry of each final state may be written at first order as [17] 

AAcP ^ agp(7r-7r+) + ^ < ^ > ^g^^ (g) 

T 

where A < t > denotes the difference of the mean proper time of the two final states and r is the true D'^ lifetime. 
Due to the difference in proper-time acceptance between the and 7r^7r+ samples A < t > /t = O.lOiO.Ol. 

Although the measured value of AAcp includes a residual 10% of the mode independent indirect CP asymmetry, 
AAcp is primarily sensitive to direct CPV. Thus the measurement of AAcp and Ap are complementary in the 
search for CP violation. The current knowledge of these measurements leads to an agreement with the no CP 
violation hypothesis with a C.L. of 20% [18] . 



5. Conclusion 

The first measurements at LHCb for search for CP violation in the charm sector are competitive with the 
results of the B-factories, even though only a total integrated luminosity of 37 pb~^ collected in 2010 is used. 
The search for direct CP violation in £)+ — >■ K~K~^Tr^ decays with a method based on the study of Dalitz plots 
indicate no evidence of CP violation. The measurement of indirect CP violation is performed on singly Cabibbo 
suppressed two-body charm decays through the study of the asymmetry of the proper-time. It is evaluated to 
be Ar = (-5.9 ± 5.9stat ± 2.1syst)- The time-integrated CP asymmetries of D° K~K+ and £>° 7r~7r+ 
decays is measured to be AAcp = (—0.28 ±0.70 ±0.25)%. In addition to these measurements, many others are 
under way, e.g. in 2-body decays the measurement of the mixing parameters using doubly Cabibbo-suppressed 
— >■ K^n~ decays. Significant improvements in the precision is expected with the large data set collected in 
2011 with an expected integrated luminosity of about 1 fb^^. 
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